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ABSTRACT 

Context. The acoustic resonator is an important model for explaining the three-minute oscillations in the chromosphere above sunspot 
umbrae. The steep temperature gradients at the photosphere and transition region provide the cavity for the acoustic resonator, which 
allows waves to be both partially transmitted and partially reflected. 

Aims. In this paper, a new method of estimating the size and temperature profile of the chromospheric cavity above a sunspot umbra 
is developed. 

Methods. The magnetic field above umbrae is modelled numerically in 1.5D with slow magnetoacoustic wave trains travelling along 
magnetic fieldlines. Resonances are driven by applying the random noise of three different colours—white, pink and brown—as small 
velocity perturbations to the upper convection zone. Energy escapes the resonating cavity and generates wave trains moving into the 
corona. Line of sight (LOS) integration is also performed to determine the observable spectra through SDO/AIA. 

Results. The numerical results show that the gradient of the coronal spectra is directly correlated with the chromosperic temperature 
configuration. As the chromospheric cavity size increases, the spectral gradient becomes shallower. When LOS integrations is per¬ 
formed, the resulting spectra demonstrate a broadband of excited frequencies that is correlated with the chromospheric cavity size. 
The broadband of excited frequencies becomes narrower as the chromospheric cavity size increases. 

Conclusions. These two results provide a potentially useful diagnostic for the chromospheric temperature profile by considering 
coronal velocity oscillations. 
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1. Introduction 


There are three main categories for the oscillations that occur 
in and around sunspots: five-minute oscillations at photospheric 
levels, three-minute oscillations above the umbra, and running 
: away from the urnbra a long penumbral structures 


waves moving a 
iBogdanl2000t[l 


_ : iBogdan & Judgel200^. Many observations s how 

the existence of three-minute waves. liVtarsh & Walshl (l2006l) ob¬ 
served these oscillations through both the 171 A TRACE and 
the SOHO/CDS instruments. The three-minute oscillations have 
also been observed at the transition region in the microwave 
band as a modulation of gyroreson ant emission (Shibasakil 
20011) and in u ltravio let wavelengths (iDe Moortel et^.ll2006h . 
Thomas et al.l (1 19871) observed multiple peaks in the three- 
minute band (4.5-10 mHz) from both space and ground in- 
str uments. A simil a r exc ited range of frequencies is present 
in iReznikova et al.l (|2012|) . A recent overview of magnetohy- 
drodynamic waves and coronal seismology can be found in 
De Moortel & Nakariakov|([2012h . 


There are three competing mechanisms to explain the source 
of the three-minute oscillations: excitation of the cut-off fre¬ 
quency, the wake of propagating shock wa ves, and a chromo¬ 
spheric resonator. lEleck & Schmitz! (1199 ll) suggested that the 
three-minute oscillations are due to a basic physical effect: 
the excitation of waves at the cut-off frequency. They numeri¬ 
cally modelled a piston driving waves at the photosphere and 
a simplified isothermal atmosphere that ignored many com¬ 
plexities including temperature gradients an d non-linearities. 
The m odel produced three-minute oscillations. iFleck & Schmit3 
(Il99ll) also investigated a VALC profile with similar results. 


Their explanation was that waves above the cut-off frequency 
are free to propagate upwards, increasing in amplitude, whereas 
waves below the cut-off frequency are damped. This is however, 
insufficient to expl ain the observed ampl ification of the three- 
minute oscillations (IStangalini et al.ll2012l) . 

It has also been suggested that the three-minute o scillations 
arise from the wake of propagating shock waves, (see lHansteenI 
1 19971) . The se shocks are als o prese nt in the ID numerical sim¬ 
ulations of iBard & CarlsonI (l2010l) . However, ID simulations 
of propagating shock wav es leads to unrealistic shock merging 
(lUlmschneider et^l2005h . 

An alternate explanation for the existence of three-minute 
oscillations above sunspot umbrae is the presence of an acoustic 
chromospheric resonator. Above the chromosphere in the tran¬ 
sition region and below it in the photosphere, there are large 
temperature gradients that provide the semi-permeable bound¬ 
aries necessary for resonances o f slow magnetoac oustic waves to 
occur inside the chromosphere (IZhugzhdall2008l) . When a wave 
encounters these large temperature gradients, part of the wave 
is transmitted and part is reflected. We n ote that the model set¬ 
up in the previously m entioned models dBard & Carlsonll2010t 
iFleck & Schmitzl[l99ll) prohibits any chromospheric resonator 
since the waves are driven using a piston at the photosphere. 
IZhugzhd^ (I 2 OO 8 I) also investigated how frequencies above and 
below the cut-off frequency behave in this resonator by using a 
temperature that changes non-monotonically with height. Waves 
above the cut-off frequency are partially reflected due to the 
sharp temperature gradients at the photosphere and transition re¬ 
gion in the sunspot atmosphere. Waves below the acoustic cut¬ 
off frequency cause the chromosphere to resonate between the 
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transition region and the point at which the cut-off frequ ency 
drops below the wave frequency (iTarovan & Erdelvll2008h . The 
resonances cause oscillations that leak energy through the tran¬ 
sition region int o the corona . This c auses regular waves to travel 
into the corona. iBotha et ^ (1201 Ih performed numerical simu¬ 
lations of a chromospheric resonator with a perturbation in the 
form of a single pulse. Two such pulses were considered, two 
and five minutes, propagating along a magnetic field line above 
a sunspot umbra. The results numerically proved the existence 
of an acoustic resonator and the production of the familiar three- 
minute oscillations. 

Observational results show that the amplitude of the three 
minute oscillation is not e xplained well enough b y the excitation 
of the cut-off frequency (IStangalini et al.ll2012h . F urthermore, 
many observations show a harmonic wave str ucture (lYuan et al.l 
1201 ll; iReznikova et al.ll2012l: ISvch et aTll201 Ih as opposed to the 
saw-tooth signal required for the shock wave model. As such, 
the chromospheric resonator model is a promising explanation 
for the three-minute oscillations. 

In this paper, perturbation noise is applied in the upper con¬ 
vection zone and allowed to propagate along the magnetic field 
line above a sunspot umbra. The wave trains will be interpreted 
as slow magnetoacous tic waves propagating with the local sound 
speed (lRobertsll2006l) . The applied noise is of the form 1 / f^, 
corre sponding to solar granulation noise (iRabello-Soares et al.l 
Il997h . 

This paper investigates the effect of different stochastic noise 
signatures on the chromospheric resonator. The effect of differ¬ 
ent chromospheric temperature conhgurations on the resulting 
coronal spectra is considered (Section |2l. The findings indicate 
a novel method for estimating the chromospheric temperature 
configuration based upon the frequency spectra present in the 
wave trains that escape from the resonating chromosphere into 
the corona. Line-of-sight (LOS) integration is then performed 
on the data using several SDO/AIA response functions (Section 

HJ. 


2. Model 

A 1.5D numerical investigation was performed to simulate the 
effects of different noise signatures and chromospheric temper¬ 
ature configurations on the coronal frequency spectra. This in¬ 
volves the generation of noise, which is described in Section 
12.31 Summaries of the numerical configuration and the post¬ 
processing method are also provided. 


2.1. Numerical set-up 

The ideal nonlinear compressible magnetohydrodynamic 
(MHD) equation s are implemented n umerically using Lare2D, 
as presented by lArber et aP (l2001ll . These equations are as 
follows: 


up 

dt 
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where p, P, T, and v are the plasma density, pressure, temper¬ 
ature, and velocity, respectively. Here, J is current density, B is 
magnetic field, E is electric field, and g is gravity. The internal 
energy e is given by Equation (0. Here, - l.^rrip where nip 
is the mass of a proton, ks is the Boltzmann constant, pQ is the 
permeability of free space, and y - 5/3 is the ratio of specific 
heats. 

The normalisation values for length, magnetic field, and den¬ 
sity are Lq = 150 x 10^ m, Bq — 0.12 T, and po = 2.7 x 10“"* 
kg m^-*. Gravity is constant as 274 m s~^ for the majority of 
the domain, however above approximately 44 Mm the gravity is 
gradually reduced to zero on the upper boundary. This is to allow 
for the specification of an open boundary and prevent reflections 
from the upper boundary. 

The velocities along a unidirectional magnetic field line 
above the umbra, perpendicular to the photospheric plane, are 
considered. The magnetic field acts as a waveguide for slow 
magnetoacoustic waves propagating into the corona. These slow 
waves propagate parallel to the magnetic field at the acoustic 
speed, which is independent of magnetic field strength. This al¬ 
lows us to choose in our 1.5D model a constant magnetic field 
strength of 10 Gauss, in line with the chromospheric strength of 
a sunspot umbral magnetic field. 

The full equations are solved in the three Cartesian coordi¬ 
nates {x, y, z) with y taken as the vertical direction. The numerical 
grid contains 1024 cells in this direction. The two horizontal di¬ 
rections {x,z} are invariant, making the simulation 1.5D. There 
are four cells in the invariant x direction. A convergence test 
of doubling the grid resolution in y resulted in no new features 
hence 1024 cells are deemed sufficient. The grid spans the range 
-68 Mm to 68 Mm, leading to a vertical cell size of 132,927 m. 
The chromosphere is positioned in the middle of the computa¬ 
tional domain. The large domain size is used to prevent bound¬ 
aries from interfering with the resonator. The end time is set as 
10"* seconds. 

The numerical dissipation in the cod e was tested by la unch- 
ing a slow pulse in a uniform plasma (e.g. lYuan et al.l20T^ using 
a range of values: p = 10“^, 10“**, and 10“*^ kg m“^ and 7" = 0.4 
MK, and 1 MK. The maximum amplitude and full width at half 
maximum were consistent to four and seven significant figures, 
respectively when calculated at the start and end of the domain. 
Therefore the numerical algorithm used here is not significantly 
dispersive. 
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Fig. 1: Initial 
profiles used 
height 0. 



temperature (solid line) and density (dashed line) 
in the simulation. The photosphere is located at 


Time, s 

Fig. 3: Example time series of the different noise sources, which 
are applied to the lower boundary as slow magnetoacoustic ve¬ 
locity perturbations. 


2.2. Temperature profiles 

The coronal temperatu re profile is obtained from the model by 
lAvrett & LoesejfcoOSi) . with the chromospheric temperature be¬ 
ing replaced by the investigated profiles. To prevent velocity per¬ 
turbations being reflected from the upper boundary, it is neces¬ 
sary that the temperature gradient at the top of the computational 
domain is zero. Therefore, from approximately 44 Mm above the 
photosphere and higher the temperature is flattened. This has no 
effect on the simulation results and the data are investigated at 
points below this region. The temperature profile is shown in 
Figure [1] 

For the convection zone, a polytrope temperature profile is 
used. The lower velocity boundary is specified directly from the 
input noise seed. By setting the lower boundary far away from 
the resonator, downwards propagating waves are damped and do 
not reflect from the lower boundary. 

There are several theoretical models availabl e for the struc¬ 
ture of the chromosphere above a sunsp ot dStaudel Il98lt 
iLites & Skumanichlll982t iMaltbv et al.lll98^ as shown in Fig¬ 
ure |2^ These theoretical models suggest that the chromospheric 
cavity size is between 1.5 and 1.8 Mm. Here two profiles are 
investigated, a flat and a step profile (Figure l2Ft and they have 
a chromospheric cavity size a + b.\r\ the numerical model, a is 
the size of the lower chromospheric temperature plateau and b is 
the size of the upper chromospheric temperature plateau. The flat 
profile allows for consideration of an idealised scenario, whereas 
the step profile is closer to the theoretical models with a lower 
and an upper chromospheric temperature plateau, as well as a 
temperature jump at the mid-chromosphere. 

At both the photosphere and transition regions there is a steep 
temperature gradient, as can be seen in Figures [T] |2^ and 
These act as semi-permeable boundaries that allow energy to be 
partially reflected and partially transmitted. The energy partially 
reflected back into the chromosphere results in resonance in the 
chromospheric c avity that exists between the photo sphere and 
transition region (IZhugzhdall2008l : [Botha et aTIliol ll) . 

In iBotha et al.l (1201 Ih . where a single pulse was used to en¬ 
ergise the chromospheric cavity, the resonator found its natural 
frequencies. In the present study, a continuous noisy source of 


energy enters the chromospheric cavity. This results in a broad 
energisation of a range of frequencies, as will be discussed in 
Section 13.31 In addition, energisation is caused by the partial 
reflection and partial transmission of the energy content at the 
photosphere and at the transition region. 

Multiple simulations were performed, varying the chromo¬ 
spheric cavity size and velocity source. This changes the energy 
content and energy distribution of the oscillations present in the 
resonator. 

2.3. Velocity initialisation 

The resonator will be perturbed by applying random noise in the 
upper convection zone, 68 Mm below the photosphere. Granula¬ 
tion noise follows a 1 //^ distribution wher e the power is propor¬ 
tional to the reciprocal of the frequency (iRabello-Soares et~an 
Il997h . Three types of noise will be considered; white (uniform 
distribution), pink (1//) and brown (1//^). This allows inves¬ 
tigation into the effect of different perturbations. The velocity 
signal applied at the source is therefore of the form 

oo 

v(f) = ^ an sin(n7rf -H p„) (9) 

n=0 

where the constants a„ is determined by the noise colour and /3„ 
is a random phase shift. 

White noise is the standard random noise where there is 
equal power in all frequencies. However since both pink and 
brown noise decline in power as the frequencies increase, these 
two are biased towards lower frequencies. Examples of the three 
types of noise are shown in Eigure|3] 

The noise samples will be generated using a Alter method. 
This involves generating a white noise sample. The white noise 
sample is then converted to spectral space via a fast Eourier 
transform (EET). The spectra are multiplied by the desired fre¬ 
quency distribution and then the signal is converted ba ck to phys¬ 
ical sp ace with an inyerse FET. Radial basis functions dBuhmannl 
l2000l) are then used to smooth the signal. An FFT is then per¬ 
formed again to this smoothed signal and the gradient is calcu¬ 
lated to ensure that it is within acceptable bounds (+0.1) of the 
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Height, Mm Height, Mm 


(a) 


(b) 


Fig. 2: Chrom ospheric temperature profiles, (a) Diffe rent theoretical temperatur e profiles for the chro mosphere abo ve a sunspot; 
solid line from iMaltbv et al.l (Il986l) . dashed line from iLites & Skumanichl(ll982h . and dotted line from lStaudel(ll98lh . (b) Flat and 
step chromospheric temperature profiles with the photosphere at 0 and the transition region at 1.8 Mm. This chromospheric profile 
is shown in relation to the photospheric and coronal temperature profiles in Figure [T] The step profile is a simplification of the 
theoretical temperature models shown in Figure |2a] 



Fig. 4; Frequency spectrum of the white noise sample shown in 
Figure |3] The pink and brown noise show similar spectra, with 
gradients of -1 and -2 respectively. 


desired gradient. Due to the stochastic nature of noise, multiple 
samples of each colour were tested to ensure reliable conclusions 
can be deduced. To have a fully resolved spectra would require 
infinite random values to be generated at the initial stage. How¬ 
ever, this is not possible. Here 1024 random numbers are used 
and the signal is well resolved up to approximately 50 mHz. 
A signal resolved to a higher frequency, 100 mHz, was tested. 
There was no difference in the results and the lower cut-off is 
used throughout this paper, i.e. 50 mHz. An example of the fre¬ 
quency spectrum for white noise is given in Figure |4l The ve¬ 
locity is scaled to be between +1 ms“* at the lower boundary, 
68 Mm below the photosphere, allowing for the consideration of 
small perturbations. 

The applied velocity corresponds to granulation noise. How¬ 
ever the source is far below the height at which one would expect 



Fig. 5: Snapshot of the evolution of velocity across the domain at 
one time instance. The source is located at -68 Mm and the pho¬ 
tosphere is at 0 Mm. Large amplitude velocity is present in the 
chromosphere due to the high energy content of the resonating 
cavity. Wave trains propagating into the corona are reasonably 
well structured. Time evolution of velocity at different points is 
shown in Figure|6] For this plot a flat chromospheric temperature 
profile of depth 1.8 Mm was used with a white noise velocity 
source. 


granules to form. The justification for this is that the nature of the 
upwards propagating velocity perturbation from the source does 
not change in the model’s convection zone, i.e. the distribution 
of noise being fed into the chromospheric resonator is the same 
as the distribution at the source, only the amplitude changes. The 
lower boundary is set far from the photosphere (-68 Mm) to al¬ 
low the downwards propagating waves from the resonator to dis¬ 
perse naturally without reflection from the lower boundary. 
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Fig. 6: Plots of velocity/sound speed vs time for different heights 
H above the photosphere, (a) and (b) are in the corona, (c) is lo¬ 
cated in the upper chromosphere, (d) in the lower chromosphere. 
The model set up for this plot was a flat chromospheric tempera¬ 
ture profile of depth 1.8 Mm with a white noise velocity source. 


3. Transmission into the corona 

3.1. Velocity propagation from the chromosphere to the 
corona 

The behaviour of the velocity perturbations at different points 
in the solar atmosphere is shown in Figures |5] |6] and |7] Figure 
|5] shows a snapshot of the velocity across the domain. Below 
the photosphere the velocity is too small to be seen on the plot. 
However there are wave trains propagating both towards the res¬ 
onator, from the velocity source applied at the lower boundary, 
and away from the resonator towards the lower boundary, due to 
downwards partially transmitted velocity perturbations. 

A time series of the velocity perturbations measured at a 
point in the lower and upper chromosphere is shown in Figures 
|6]and|2] The wave trains propagating into the corona appear well 
structured, i.e. well resolved with no shocks forming, with sev¬ 
eral higher frequencies present. This is reflected in the frequency 
spectra. 

Figure |6] shows the time evolution of the velocity at different 
spatial points. The velocity amplitude, with respect to the lo¬ 
cal sound speed, increases throughout the chromosphere. When 
energy is transmitted through the temperature boundary at the 
transition region, there is a decrease in power since the energy 
is partially reflected back into the chromosphere. Again the sig¬ 
nal appears to be mostly regular with several higher frequencies 
present, identifiable from the small scale fluctuations. Note that 
the two coronal figures, (a) and (b) in Figure |6l show the same 
velocity perturbation shifted by time. This shows that once the 
signal has entered the corona, there are no significant changes 
to the velocity perturbation. There is no evidence of any shocks 
occuri'ing and the signal is well resolved. The frequency regimes 
discussed in Section [33] are all present when sampled at various 
points in and above the transition region. 



Time, s 


Fig. 7; Plots of velocity/sound speed vs time for different heights 
H above the photosphere, (a), (b), (c) and (d) are at the same 
locations as in Figure |6l A step temperature profile was used 
with a - b - 0.9 Mm for both the upper and lower plateaux. 
The perturbation noise was white. 


3.2. Frequency analysis 

An FFT was performed to analyse the frequencies present in the 
corona for each simulation. This involved considering the time 
series velocity fluctuations through a single point in space. A line 
has then been fitted to the FFT to quantify the relation between 
power and frequency using the well established method of least 
squares. The gradient of this line is then used to consider the 
effect the chromospheric temperature configuration has on the 
signal propagating through the corona. Velocity data are sampled 
at intervals of 0.25 seconds, leading to a Nyquist frequency of 
2000 mHz. 


3.3. Coronal frequency regimes 

The frequency for each simulation is sampled through time in 
the corona at a point 25 Mm above the start of the transition 
region. Four main regions are observed in the spectra for all the 
simulations, as shown in FiguresOand|9]for the flat and step pro¬ 
files respectively; all except the lower cut-off frequency are due 
to the resonating behaviour. There is a general trend of increased 
spectral power for larger chromospheres. This is expected since 
there is a larger region for the resonances to oscillate in, which 
provides increased amplitude and hence higher power. Note also 
that there is a level of uncertainty present in the data due to the 
fact that the system is being driven by stochastic noise. This 
arises because the noise is generated using a finite number of 
samples resulting in slight deviations from the expected power at 
each frequency. At double the grid resolution, the four regimes 
were still present at the same frequencies. 

Region I: Lower Cut-off Below approximately 4 mHz there 
is a very low power region (Figure [8] region I). This is due to 
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Fig. 8: Four frequency regimes present in the spectra. I is the 
lower cut-off due to the acoustic cut-off frequency (below 4 
mFIz). 11 is the broad peak (4-17 mFIz). Ill is the gradient de¬ 
cline (17-200 mHz). IV is the upper cut-off (above 200 mHz). 
This spectrum was generated for a flat chromospheric profile of 
depth 1.8 Mm with a white noise velocity source. The line in¬ 
dicates the gradient of the power in region III, which is -1.15 in 
this case. 



Fig. 9: Four frequency regimes present in the spectra. Regions 
I, II, III and IV are the same as in Figure 0 generated for a 
step chromospheric profile with upper and lower plateau depths 
of a = b = 0.9 Mm and a white noise velocity source. AT is 
2600 K. The line indicates the gradient of the power in region 
III, which is -2.62 in this case. 


the acoustic cut-off frequency ojc that allows waves to propagate 
upwards with frequencies 


to > tOc 


Cs 

2H 



( 10 ) 


where H = y is the ra tio of specific heats, Cj is the 

sound speed, and z is the height (lRobertsll2004t) . Typical values 
of the acoustic cut-off frequency are in the range 3.5 - 5.2 mHz 


(e.g.,lBel & Lerovlll977l: iMcIntosh & Jefferi^l2006l: lYuan et all 

I2014h . 

Region II: Broad Peak Between approximately 4 and 17 mHz 
the power is relatively flat (Figure^ region II). This corresponds 
to a broad band of frequencies with oscillation times of be¬ 
tween a pproximately one-m i nute to four-minutes. Following the 
work of iBeckers & Schultzl dl972h . several papers have inves¬ 
tigated t he three-minute oscillat ions above sunspot umbra. For 
example. [Reznikova et al.l (12012h present spectra for oscillations 
above sunspot umbrae through various filters which show an 
excited range of fre quencies between approximately 6-8 mHz. 
iThomas et al.l (Il987h makes reference to a three-minute band of 
frequencies between 4.5 and 10 mHz containing multiple peaks. 

This broad peak is present in the same frequency range re¬ 
gardless of temperature configuration or perturbation noise. No 
discrete harmonics are visible in the spectra. This does not mean 
that the harmonics are not present; they are not visible due to the 
power distribution in the frequency range where the harmonics 
are situated. 

Region III: Gradient Deciine The next region (Figure 0 re¬ 
gion III) between 17 and 200 mHz, is a region where the power 
is proportional to 1 //“, which corresponds to a linear decline of 
gradient a when plotted on log-log axes. This region is found 
to vary in gradient when the temperature profile is modified and 
provides the diagnostic tool for the main conclusions drawn in 
this paper. This is discussed in more detail in the Sections (3.4.II 
and l3.4.2l Confidence intervals were calculated to ensure the fit¬ 
ted line was representative of the data. The standard deviation 
was small in all cases, ranging from 0.041 to 0.052, resulting in 
very narrow confidence intervals, identical to 3 decimal places. 
For one such line, the 95% confidence interval is -0.9649 to 
-0.9646, with the calculated value being -0.9647. 

Region IV: Upper Cut-off Above approximately 200 mHz 
(Figure [8] region IV) the power drops significantly. The cut¬ 
off prevailed at the same frequency with a higher resolved in¬ 
put signal indicating that this upper cut-off is due to neither 
the input signal (Section 12.3b nor the sampling rate (Section 
Ira . This is a fluid model which does not account for the ion- 
gyrofrequency. However the ion -gyration fre quency is approxi¬ 
mately Q, = 9.6 X 10"^ rad s“^ (lPriestll20T4l) . This is far larger 
than the frequencies simulated in this model, therefore the ob¬ 
tained frequencies are not beyond the physical limitations of the 
fluid model. This suggests a physical phenomena for the upper 
cut-off unidentified at present. Due to the high frequency where 
the cut-off occurs, it cannot be resolyed in obseryations at the 
present time. 

3.4. Varying chromospheric cavity configuration 
3.4.1. Flat Profile 

Initially the temperature configuration of the chromosphere was 
assumed to be uniform, as shown in Figure |2b] with AT = 0. 
This was to inyestigate the behayiour in an idealised scenario. 
Here the parameter space consists of the chromospheric cay- 
ity size {a + b in Figure l2bT l and the perturbation noise colour. 
Three noise samples for each noise colour were used and a range 
of chromospheric cayity sizes between 0.66 and 2.13 Mm were 
tested. The gradient of the aforementioned frequency range (Fig- 
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Chromospheric cavity size, Mm 



(a) 


(b) 


Fig. 10: Gradients vs chromospheric cavity size for three samples of each noise colour (a), and averaged across each noise colour 
(b). The gradients were obtained in region III of Figure [8] 


ure 0 region III) was calculated for each simulation. Note that 
due to the stochastic nature of noise, different samples of the 
same noise colour create slightly different results, indicated by 
the variations present in Figure fTOal Multiple seeds were tested 
for this reason. 

For clarity, the gradients were averaged for each noise colour 
and are shown in Figure II Obi This effectively removes some of 
the uncertainty present when considering noise. The full set of 
results is shown in Figure llOal showing the extent of random¬ 
ness. The variation in the gradient is approximately +0.5 for nar¬ 
row chromospheres of cavity size less than 1.3 Mm, whereas for 
large chromospheres of cavity size larger than 1.8 Mm there is 
virtually no variation in gradient across samples. 

For narrow chromospheric cavities, less than 1.3 Mm, all 
three noise colours appear to have the same output gradient in 
their spectra, as shown in Figures [TOa] and [TOb] There is a sharp 
change in output gradient for chromospheric cavity sizes starting 
from around 1.3-1.4 Mm and ending at approximately 1.8 Mm. 
Here the gradient rapidly goes from steep, a » -3, to relatively 
flat a -1. Since this large change in gradient is present, it 
may be possible to estimate the depth of the chromosphere from 
the spectra sampled in the corona at 25 Mm above the transition 
region. 

The dependence of the power gradient on the chromospheric 
cavity size can be explained by considering the energy content of 
the chromosphere. As the chromospheric cavity size increases, 
the total energy content of the resonating waves increases. This 
allows for energy to be more evenly distributed among the higher 
frequencies when partially transmitted into the corona. For a 
thinner chromosphere, there is a low total energy content and as 
such only the lower frequencies become energised. As the chro¬ 
mosphere cavity increases in size, the energy content increases 
and there is more energy in the higher frequencies, relative to 
the lower frequencies, i.e. a decrease in the slope of the power 
gradient in the FFT spectra (Figure fTObl i. 

With regards to noise, there does not appear to be a sig¬ 
nificant difference between noise colours. Whilst the location 
of the sharp gradient change in the spectra measured for nar¬ 
row chromospheres appears slightly different in Figure [TOb] for 
different noise colours, there is no distinct relationship and the 
variation in the location of the sharp change is small, occuri'ing 



0.0 . 

0.0 0.5 1.0 1.5 2.0 

Upper Plate, Mm 


Fig. 11: Contours of the gradient in the obtained coronal spectra 
(region III of Figure [8]l generated by white noise and the step 
temperature profile. The lower plateau sizes are given by values 
of a and the upper plateau sizes by b in Figure The black 
marker on the vertical axis indicates the steep gradient change 
for the flat profile, as seen in Figure fTObl 


between approximately 1.3 and 1.4 Mm of chromospheric thick¬ 
ness. Therefore the discrepancy is most likely due to the stochas- 
ticity present. 

For large chromosphere depths, above 1.8 Mm, the output 
gradient becomes fairly constant for all three noise colours, as 
shown in Figures fTO^ and fTObl 

3.4.2. Step profile 

Next, the chromospheric temperature configuration is replaced 
with a step profile. This consists of two uniform temperature 
regions or plateaux, a and b, separated by a small temperature 
jump, AT as shown in Figure |23 This is similar to theoreti¬ 
cal chromospheric temperature models (Figure |2a|l and provides 
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a more realistic configuration. Thus the parameter space here 
consists of 4 variables; lower plateau size a, upper plateau size 
b, temperature jump AT, and noise colours. In this section, the 
temperature jump will be fixed at 2600 K since this agrees with 
the theoretical models (Figure|2a]i. Varying the temperature jump 
will be investigated in Section l33] 

The depth of both the lower (a) and upper (b) plateau varied 
in the range 0.13 Mm < a,b < 2.0 Mm, under the constraint that 
the total depth 0.8 Mm < a + b < 2.13 Mm. This provides a 
maximum total chromospheric cavity size of approximately the 
largest value from the theoretical models, i.e. 2 Mm as shown in 
Figure |2a] The minimum chromospheric cavity size is 0.8 Mm. 
The three types of noise were used as in Section [3 .4. II however 
on this occasion only one sample of each colour was used. The 
gradient of the coronal spectra was calculated as before in region 
III of Figure 0 The results can be seen in Figure [TT] for white 
noise. 

A similar trend is present here as with the flat profile: there 
is a steep change in spectral gradient between two regions of 
near-constant spectral gradient, as shown in Figure fTObl The lo¬ 
cation of the steep gradient change in Figure [TT] is independent 
of the noise source. An interesting feature is that the location of 
the steep spectral gradient change varies with plateau sizes. The 
location of the sharp change approximately follows the equa¬ 
tion a = -Q.5b H- 1.5 for lower plateau size a and upper plateau 
size b. As the size of the upper plateau (b) increases, the size 
of the lower plateau (a) required for the sharp gradient change 
decreases but at a slower rate. This implies that the dimension 
of the lower plateau a is more important in producing this steep 
change than the dimension of the upper plateau b. The marker 
on the vertical axis in this plot indicates the location of the steep 
change for the flat profile. 

For a small lower plateau , there is near constant spectral 
gradient, similar to the results from the flat profile. Figure llOal 
In Figure [TT] there is some variation present which is due to the 
stochastic element of these simulations. One would expect that 
averaging over several samples of the same noise colour would 
result in a smoother result, as with the flat profile. Figure [TOb] 

As in the case of the flat profile, there is very little difference 
in the location of the sharp gradient change between the results 
for different noise colours. Any discrepancies are of the same 
order as those found for the flat profile (Section [3.4.1b . For this 
reason, only the white noise results are presented. It can be seen 
that the trend for the flat profile (Figure [TObT i continues the trend 
present for the step profile. Figure [TT] 

3.5. Varying the mid-chromospheric temperature jump 

A set of simulations was performed to investigate the importance 
of the temperature jump AT between the two uniform regions. 
For this, only a few chromospheric cavity sizes were tested, 0.6 
Mm < a-\-b < 1.2 Mm, and three temperature jumps were con¬ 
sidered; AT - 1600, 2600 and 3600 K. The results are shown 
in Figure [12] There appears to be no significant difference be¬ 
tween the temperature jumps, indicating that the variable AT is 
less important than the total chromospheric cavity size a b. 


4. Chromospheric seismology 

The steep change in spectral gradient can be clearly seen in Fig¬ 
ures [8] and [9] Regions I and II in these figures are very similar. 
The diagnostic in this paper is demonstrated by considering re¬ 
gion III of these two plots. In Figure [8] the chromospheric tem- 



Fig. 12: Gradients of the spectra (Figure^ region III) obtained 
with different temperature jumps at the mid chromosphere, i.e. 
with different values of AT as shown in Figure [2b] 1600K, 
2600K, 3600K. The chromospheric cavity size a b increases 
with a - b. 

perature profile leads to a shallow gradient, whereas the temper¬ 
ature configuration in Figure 0 leads to a steep gradient. One 
can clearly see the stark differences between these two plots for 
region III. 

The results obtained give rise to a new potentially useful di¬ 
agnostic for chromospheric temperature configurations. Given 
the velocity spectra in the corona, it is possible to give upper and 
lower limits to the chromospheric cavity size. This diagnostic 
arises from the presence of the steep change in gradient present 
in all cases considered here, visible in Figures II Obi and [TT] For 
example, consider a case where the spectral gradient is steep. 
From Figure [TT]it can be seen that the size of the chromospheric 
cavity is bounded above by the steep gradient change, which fol¬ 
lows the line a = 0.55 H- 1.5. This implies limits on the size of 
the lower plateau 0 < a < 1.5 Mm and upper plateau 0 < 5 < 3 
Mm, subject to the constraint that a ■+■ 0.55 <1.5. 

5. Spectral shape and its height above the 
photosphere 

In addition to the spectra in the corona, the photospheric spectra 
were also analysed for white noise. There is a reasonable flat 
spectra followed by a sudden drop at approximately 20 mHz, 
as shown in Figure [13] Again this sudden drop in power was 
found to be independent of the input wave resolution. A possible 
explanation for this is that the downwards propagating waves 
due to the resonator dwarf the driven upward waves in power at 
this point. Note that this does not mean the upward propagating 
waves are not present, merely that they are not obvious in the 
FFT since the waves propagating from the resonator are far more 
powerful. 

It was also found that the spectra vary through the chromo¬ 
sphere. In particular, either side of a temperature jump i.e. the 
mid-chromospheric jump (AT in FigurelZbl) and at the transition 
region, the higher frequencies in the spectra change behaviour. 
The upper spectral cut-off in Figure [T3] is very steep. However 
if a mid-chromosphere temperature jump is present, the signal 
in the upper chromosphere has a more gradual cut-off. This de- 
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Fig. 15. Line of sight integration of the coro¬ 
nal velocity perturbations for 171 A for four 
different sizes of the chromospheric cavity. The 
power spectra are normalised to their respec¬ 
tive maximum values for chromospheric cavity 
sizes (a-rb in Figure l2bl l of the temperature pro¬ 
file given by 0.67, 1.20, 1.60 and 2.13 Mm re¬ 
spectively. 



Fig. 13: Frequency spectrum sampled at the photosphere, i.e. at 
height 0. 


velops into the gradient and cut-off regions III and IV seen in 
Figures [8] and 0 above the transition region. Also, as expected, 
the acoustic cut-off frequency in Region I becomes more promi¬ 
nent towards the transition region. 


94A 


131A 



Height, Mm Height, Mm 

17iA 



Height, Mm 


Fig. 14: Normalised response function RiT) plotted against 
height using three SDO/AIA filters at time t - 0. This is recal¬ 
culated at each time iteration for the LOS integration to account 
for temperature variations. The corresponding initial tempera¬ 
ture profile is shown in Figure[T] 


6. Limitations of the observabie coronai spectra 

In order to investigate how the findings presented in this paper 
would be observed, LOS integration is performed. The line of 
sight integration is calculated using 

I^jR(T)nldy (11) 

where I is the resultant line intensity, is the number density, 
R{T) is the response function for temperature T, and y is the 


height. Note that the mass density p is proportional to the num¬ 
ber density rig. This integration was performed from the transi¬ 
tion region to the corona. This was performed numerically using 
the trapezium rule. The temporal resolution here is one sample 
every two seconds. This corresponds to a Nyquist frequency of 
250 mHz. The response function is recalculated at each time step 
to account for temperature fluctuations. The integration was per¬ 
formed with a lower bound as the start of the transition region to 
remove the chromospheric foot-point. Only the flat profile driven 
by white noise has been considered. The line of sight integration 
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was performed for chromospheric cavity sizes between 0.66 and 
2.13 Mm. 

6.1. SDO/AIA response functions 

The initial response function R(T) as a function of height is 
shown in Figure [14] for three SDO/AIA filters (94, 131, and 
171 A). The SDO/AIA response functions were g enerated using 
the m ost recent version V6 available in ssw-idl (iBoerner et al.l 
I 2 OI 2 I) . LOS integration was performed on the three samples of 
white noise, varying the chromospheric cavity size between 5 
and 2.13 Mm for the three SDO/AIA filters shown in Figure [T4l 
The three SDO/AIA filters give similar results, hence we only 
discuss the behaviour for the 171 A channel. 

Figure [15] shows the resultant frequency spectra of the line 
intensity using the 171 A filter for chromospheric cavity sizes 
of 0.66, 1.20, 1.60 and 2.13 Mm. The spectra demonstrate the 
acoustic cut-off for frequencies less than 4 mHz, similar to the 
spectra obtained from a point source (Figures [8] and ]9] region I). 
As the chromospheric cavity size increases, the power of the un¬ 
derlying noise becomes less dominant and the highest frequency 
in the excited range decreases. In other words, the bandwidth 
containing the excited frequencies becomes narrower. For exam¬ 
ple, in Figure [ 15 ] the excited range of frequencies for a chromo¬ 
spheric cavity size of 1.595 Mm is between approximately 4 and 
18 mHz whereas for a chromospheric cavity size of 2.127 Mm 
the range is between 4 and 13 mHz. 

A line of the form 1 //“ can be fitted to the LOS spectra in 
the range 20-50 mHz, with a representing the log-log gradient. 
For the single point analysis there was a trend present where the 
gradient depended on the chromospheric cavity size. However 
only a weak trend is present in the gradients for the 94 and 171 
A filters. No discernible trend was present for the 131 A fil¬ 
ter. Variations in a between different samples of white noise are 
large with respect to the total change in a. Therefore it would be 
difficult to accurately determine the chromospheric cavity size 
given a. 

6.2. LOS integration \/s singie point anaiysis 

The LOS spectra are different to the frequency spectra in the 
single point analysis. The primary difference is the range of fre¬ 
quencies present in the broad peak. The lower limit of this region 
is the same, indicating that the acoustic cut-off frequency is de¬ 
tectable in both cases. However the upper limit of the broad peak 
is 10 mHz for the LOS spectra, but 17 mHz for the single-point 
spectra. 

The LOS integration gives different results compared to the 
single point analysis. Possible explanations are as follows: 

1. The integration of a sinusoidal function depends on the depth 
of the LOS integration and the periodicity of each individual 
component of the signal. As such, the power of the individ¬ 
ual frequencies is modified non-uniformly by performing the 
LOS integration. This results in a different power distribution 
compared to the single point analysis. 

2. The spectral lines/filters have a finite width covering a tem¬ 
perature range. This is in contrast to the single point analysis 
that isolates a single temperature value. The same holds for 
the mass density, p. 

3. The distribution of the frequency spectra changes with height 
in the transition region where the temperature gradients 
are large. As one moves into the corona, the variation in 


the shape of the frequency spectra decreases. Any spectral 
lines/filters including contributions from the transition re¬ 
gion will be affected. The single point analysis results dis¬ 
cussed in this paper are obtained in the corona. Therefore, 
the weighted transition region will enhance the differences 
between the single point analysis and the LOS integration. 

6.3. Angie of LOS integration 

In the 1.5D model the LOS integration is performed parallel to 
the magnetic field lines. Performing the integration perpendicu¬ 
lar to the magnetic field will yield the single point result (since 
it is invariant in this direction). Integrating at an angle produces 
the same results as integrating parallel to the magnetic field line, 
since, even though the path length appears longer, it intercepts 
the same discrete values in the same order. 

6.4. Comparison with observations 

The frequency spectra obtained in Figur e [TSl are v e ry sim ilar 
to observational papers as can be seen in iTian et alJ d2014li : an 
excited range of fr equencies is present between approximately 
4.5 and 10 mHz (iThomas et al.l [19^ iReznikova et al.l I 2 OI 2 I : 
iTian et a02014ll . For wider chromospheres, this is more distinct 
and the excited range is reasonably narrow. However for smaller 
chromospheres the LOS spectra are far more noisy and whilst 
there are a few strong peaks, the power does not decay rapidly, 
as can be seen in Figure[T5] The bandwidth of the LOS frequency 
spectra can be used to estimate the chromospheric cavity size, as 
was discussed in Section 16.11 For exam ple the 4.5 to 10 mHz 
range observed by iThomas et al.l (Il987h would indicate a chro¬ 
mospheric cavity size of greater than 2.13 Mm. 

7. Conclusions 

In this paper, a numerical investigation of oscillations above a 
sunspot umbra was performed. Multiple simulations were per¬ 
formed where both the chromospheric temperature profile and 
perturbation noise were varied. The resultant coronal velocity 
signature was then analysed. Different colours of noise were 
found to have very little effect on the output spectra. However, a 
strong trend is present where the gradient in the frequency spec¬ 
tra of the coronal velocity perturbations varies with regards to 
chromospheric cavity size: as the chromosphere cavity increases 
in size, the spectral gradient becomes shallower. 

The coronal spectra can potentially be used as a diagnostic 
for the chromospheric cavity size. The sharp change in gradient 
present in Figures [1 Obi and [m allows for limits to be estimated 
for the sizes of the lower and upper uniform temperature regions 
or plateaux (a and b in Figure l2bTl. These limits can be further 
reduced by estimating a maximum and minimum total chromo¬ 
spheric cavity size (a + b). 

The LOS spectra possess high power oscillations in the 
three-minute band, similar to those in observed spectra. It is 
found that the bandwidth of the excited frequencies become 
narrower as the chromospheric cavity size increases. This is a 
promising tool to estimate the chromospheric cavity size. 

There was no clear relation present between the LOS spec¬ 
tral gradient and the chromospheric cavity size when integrating 
along a magnetic field line. The LOS integration along a mag¬ 
netic field line masks the underlying complexity of the waves 
present in the single point analysis. This indicates that a very 
narrow response function is necessary to clearly observe a trend 
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along a magnetic field line. Alternatively, observing the sunspot 
at an angle to the photosphere would yield a different view where 
the waves are not integrated along the magnetic field line. As 
such, one would expect the results to approach the results ob¬ 
tained from the single point analysis as the angle approaches 90° 
with respect to the magnetic field lines. 

It is inherently difficult to measure the chromospheric cav¬ 
ity size directly. This work provides a solution by analysing the 
coronal waves above sunspots for different chromospheric cav¬ 
ity sizes via a numerical simulation. The single point results 
presented in this paper show that the gradient of the coronal 
frequency spectra is directly correlated with the chromospheric 
cavity size. This is equivalent to performing LOS integration 
perpendicular to the magnetic field line. When LOS integra¬ 
tion is performed along the magnetic field line, it is shown that 
the bandwidth of excited frequencies varies with chromospheric 
cavity size. Thus, a novel diagnostic was presented for indirectly 
estimating the chromospheric cavity size above sunspot umbrae 
for the Sun and other stars. 
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